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Abstract 

We discuss the structure of infrared singularities in on-shell QCD amplitudes with massive partons 
and present a general factorization formula in the limit of small parton masses. The factorization 
formula gives rise to an all-order exponentiation of both, the soft poles in dimensional regular- 
ization and the large collinear logarithms of the parton masses. Moreover, it provides a univer- 
sal relation between any on-shell amplitude with massive external partons and its corresponding 
massless amplitude. For the form factor of a heavy quark we present explicit results including the 
fixed-order expansion up to three loops in the small mass limit. For general scattering processes 
we show how our constructive method applies to the computation of all singularities as well as 
the constant (mass-independent) terms of a generic massive n-parton QCD amplitude up to the 
next-to-next-to-leading order corrections. 



1 Introduction 



Amplitudes for hard scattering processes in Quantum Chromodynamics (QCD) are of basic im- 
portance both for theory and phenomenology and predictions for these matrix elements have to 
include higher-order quantum corrections. They are mandatory for precision measurements of 
Standard Model parameters and critical to the determination of backgrounds for new physics phe- 
nomena. Many explicit computations of hard multi-parton processes do not only provide us with 
a wealth of information but have also helped significantly in understanding underlying principles 
such as factorization or the universal structure of collinear and infrared singularities. 

These singularities are particularly prominent for at least two reasons. First of all, the in- 
dependent knowledge of the universal limits when parton momenta become collinear or a gluon 
momentum tends to zero serves as a very strong check on any complete calculation. Secondly, 
the calculation of finite cross-sections in QCD beyond leading order has to combine consistently 
squared matrix elements with different numbers of partons in the final state. In any such formalism 
(see e.g. [1]) the individual contributions have to be suitably integrated over the available phase 
space and are usually infrared divergent. At next-to-leading order (NLO) in QCD the singular 
behavior of the corresponding amplitudes with both massive and massless partons in the final state 
has been extensively studied [2-4]. 

Research beyond NLO in the past years has been primarily focused on the calculation of mass- 
less amplitudes at next-to-next-to-leading order (NNLO), see for example [5-9] and numerous 
references therein. The progress at NNLO and investigations of the singular behavior of am- 
plitudes at higher loops [10] have significantly contributed to our understanding of their general 
structure from the view point of all-order resummations [11]. This in turn leads to predictions for 
the soft and collinear behavior of massless amplitudes at any order based on a small number of 
perturbatively calculable anomalous dimensions. 

For massive amplitudes however much less is known beyond NLO in QCD despite the fact 
that NNLO precision predictions with massive quarks are clearly needed in view of the data from 
present and the prospects of future high-energy colliders (see Refs. [12-14] for related progress). 
Prominent examples of such measurements are for instance the forward-backward asymmetry Afb 
for inclusive heavy quark production in e + e~ -annihilation [15, 16], and cross-sections for heavy 
flavor production and decays at the Tevatron and the LHC (see e.g. Ref. [17]). 

The aim of this article is a first systematic investigation of the structure of massive QCD am- 
plitudes in singular limits beyond NLO. To that end, we extend the studies of Refs. [10, 11] to 
partonic scattering processes including the presence of massive particles. The masses of the latter 
screen the divergences of the massless amplitudes and give rise to large logarithmically enhanced 
contributions of Sudakov type [18], which dominate the high energy behavior of the scattering 
amplitudes. It is precisely the structure of these large logarithms together with soft singularities 
appearing as poles in (d — 4) in d dimensions, that we wish to address here for a general non- 
Abelian SU(iV) -gauge theory such as QCD. Throughout the article we neglect power corrections 
in the parton masses m. 



1 



The outline of the paper is as follows. In Section [2] we recall the general framework for the 
factorization of n-parton amplitudes in QCD and discuss its modifications to incorporate massive 
partons. As a result we derive an extremely simple universal multiplicative relation between a 
massive amplitude in the small mass limit and its massless version. This is one main result of this 
paper. The corresponding multiplicative factor (which we call Z) can be linked to the QCD form 
factor of massive and massless partons. Next in Section[3l we specifically address the resummation 
and exponentiation of the QCD form factor for heavy quarks, which is our second main result. On 
this basis we provide in Section @] all resummation coefficients and new fixed-order expansions of 
the massive form factor up to three loops. For the resummation coefficients we observe striking 
relations between the massless and the massive case. In Section |4] we also present explicit results 
for the universal multiplicative factor Z up to two loops and discuss its relation to the perturbative 
fragmentation function of a heavy quark [19]. We argue that our formalism represents the proper 
generalization of Ref. [19] at the level of amplitudes. In Section [5] we demonstrate the predic- 
tive power of the factorization ansatz for QCD amplitudes with examples from 2 — >• n scattering 
processes, such as hadronic ff -production. There we discuss the complete structure of the soft 
and collinear singularities including the logarithmically enhanced terms to NNLO in perturbative 
QCD. We summarize in Section |6] and present some technical details in the Appendix lAl 



2 Factorization of QCD amplitudes 

We are interested in a general 2 — > n scattering processes of partons pi 

p: Pi(ki,mi,ci)+p2(k2,m2,c 2 ) -> p 3 (k 3 ,m 3n c 3 ) + . . . + p n+ 2(k n+ 2,m n+2 ,c n+2 ) , (1) 

where {/?,} denotes the set of partons (of specific flavors) with associated momenta {ki}, masses 
{m,} and color quantum numbers {c,}. The latter are in the range 1 . ..N 2 — 1 for particles in the 
adjoint (gluons) and 1 . . .N for particles in the fundamental representation (quarks) of a SU(iV)- 
gauge theory. 

The scattering amplitude M M for the process (OQ) is conveniently expressed in a basis of color 
tensors (cj)t c .y Following Ref. [11] we write M ' p ' as 

^ } (w.^.«.(^). e ) = ^ / [pl ({^},^,a s (A' 2 ),e) (cj) {ci} (2) 

= \M P ), 

where p is the renormalization scale, Q the hard scale of the process typically related to the center- 
of-mass energy, e.g. Q = y/s with s = (ki +k 2 ) 2 , and £ the parameter of dimensional regularization, 
d = 4 — 2e. The amplitude |fW p ) is a vector in the space of color tensors cj with summation over 
/ being understood. We consider M ^ at fixed values of the external parton momenta ki, thus 
kf = mf and especially kf = for massless partons. Any additional explicit dependence on the 
parton masses m ; in Eq. © is suppressed. 
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Let us start by recalling that on-shell amplitudes for massless partonic processes in d = 4 — 2e 
dimensions can be factorized into products of functions Jq , 5q P ' and H I p ' . These functions are 
called jet, soft and hard functions and are known to organize the contributions of various momen- 
tum regions relevant to the structure of the singularities in the scattering amplitude. Following 
Refs. [10, 11] we can write 

\M P ) = 4 P] (^« S (A' 2 ),e)4 P] ({^}'^ a s(A/ 2 ),e^p). (3) 

The short-distance dynamics of the hard scattering is described by which is infrared finite. 
Analogous to the decomposition in Eq. © |#p) is a vector in the space of color tensors cj. Coherent 
soft radiation arising from the overall color flow is summarized by Sq \ where we also use matrix 
notation suppressing the color indices. The function depends only on the external partons and 
collects all collinearly sensitive contributions. It is otherwise independent of the color flow. 

Given the factorization formula © one can then organize the singularity structure of any mass- 
less QCD amplitude. After the usual ultraviolet renormalization is performed, these singularities 
generally consist of two types, soft and collinear. Being of infrared origin, they are related to 
the emission of gluons with vanishing energy and to collinear parton radiation off massless hard 
partons, respectively. In this way all soft and collinear singularities in massless amplitudes are 
regularized and appear as explicit poles in £ as indicated in Eq. ©. Typically two powers of 1/e 
are generated per loop. 

When masses are introduced the picture described above gets modified. In QCD, which has 
only massless gauge bosons, the soft singularities remain as single poles in e while some of the 
collinear singularities are now screened by the mass m of the heavy fields. Nevertheless, in pres- 
ence of masses, we speak of quasi- collinear singularities [2] that exhibit logarithmic dependence 
on m. To be specific, in the present paper we will consider the small mass limit of massive QCD 
amplitudes M M such as in Eq. ©. Naturally, in this limit we require that all masses in the ampli- 
tude are either zero or equal to a common value m and much smaller than the characteristic hard 
scale Q of the reaction. Thus, in the limit Q 2 ^> m 2 we aim at organizing all poles in £ and all 
powers of In (ra), k > 0, (including mass independent terms) from the underlying factorization 
principles. 

From an alternative point of view however, the differences between a massless and a massive 
amplitude for a given physical process can also be thought of as a mere change in the regulariza- 
tion scheme. Here, the limit of small masses for any given amplitude may simply be seen as an 
alternative to working in J-dimensions in order to regulate the soft and/or collinear singularities. 
Of course, gauge invariance has to be retained. In this interpretation parton masses act as formal 
regulators and massive amplitudes in the limit m 2 <C Q 2 must share essential properties with the 
corresponding massless amplitudes. Such arguments have been previously used in Refs. [20-22] 
in the context of QED corrections to the Bhabha process. Within QCD with n t light quarks and 
one heavy flavor, this requires to properly account for the decoupling of the heavy quark. We 
will further elaborate on this point below, in particular on the relevant aspects of the decoupling 
theorem [23]. 
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Our goal is the generalization of the infrared factorization formula © of Refs. [10, 1 1] to the 
case of massive partons. To that end, we perform a similar factorization for the amplitude M ^ 
into products of functions J7 ' p ', S ^ and 9{ In the presence of a hard scale Q we can then write 
for the partonic process © 



\M V ) = ^W(%,^,a.(^),6)5lP]({* i },%,a.(^),6)|^> J (4) 



where all non-trivial mass dependence enters in the functions J7 ^ and S ^ and we neglect in H ^ 
power suppressed terms in the parton masses m. The jet function now summarizes all quasi- 
collinear contributions from the external partons. It is therefore of the form 

'§ ' as(A,2) ' £ ) = n/ [il (f ' as(A/2) ' £ )' (5) 

where 3 M denotes the jet function of each external parton p\. 

We stress that the above factorization formula © is designed to correctly reproduce the leading 
power in the hard scale Q. Moreover, as the similarity between Eqs. © and © suggests, the 
factorization is otherwise independent of details such as the partons in reaction © being massless 
or massive. However, Eq. © still contains ambiguities related to the separation of finite terms in 
3 [ p l , S ^ and 9{ M . It also contains ambiguities related to sub-leading soft terms in 3 M and S ^ . 
Following Ref. [11] we fix this remaining freedom completely by demanding that 

Q 2 rf„r.^„\ ( rr \{\(Q 2 m 2 i 



where the scalar function f ^ denotes the gauge invariant space-like form factor of a quark or gluon 
to be discussed in detail in Section [3] below. For the moment, suffice it to say that the function f ^ 
is associated to the vertex y*qq (or J*qq), of a photon y* with virtuality Q 2 , and q/q an external 
quark /anti-quark of mass m q . Likewise, for a colored parton in the adjoint SU(iV) -representation, 
the function f ^ is either obtained from the effective vertex §gg of a scalar Higgs and two massless 
gluons, or from the corresponding vertex with two gluinos g of mass trig. 

The motivation for the choice made in Eq. © above comes from the following consideration. 
Firstly, it reproduces the collinear dynamics as desired and, moreover, provides a specific prescrip- 
tion for the pure soft terms contained in the jet function. Secondly, it guarantees that the jet factor 
3 [ p l remains process-independent, while all process-dependent soft interference terms are entirely 
delegated to the soft function S ^ . We recall that the role of parton masses is to simplify screen 
the collinear singularities. Since the soft and hard functions 5 W and are insensitive to these 
collinear dynamics, being the same in the massless or the massive case (provided Q 2 ^> m 2 ), log- 
arithmically enhanced contributions of the type \n k (m) are contained solely within 3 In other 
words, we require the (massive) factorization formula © to be valid for any amplitude. Then it 
also holds for the form factors f ™ in Eq. © itself, since these are the simplest amplitudes to which 
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Eq. © can be applied with S I" - * 1 ' = 1 and H i il ~ > 'M = 1, and this choice for J7 M is also consistent 
with the corresponding massless case. 

We also want to comment briefly on evolution and exponentiation. In Eqs. © and © we 
have suppressed any additional scale dependence, which together with the renormalization group 
properties gives rise to evolution equations for j7q P ', 5q P ' and 5^. The solution of those 
evolution equations leads to an all-order exponentiation in terms of the corresponding anomalous 
dimensions, which is well known for massless partons, see e.g. Refs. [11, 24]. In the case of 
massive partons, the exponentiation of the jet function J7 M (the form factor J ^, respectively) is 
discussed in detail in Section [31 while we postpone the soft function S ^ and its solution as a 
path-ordered exponential until Section [51 

Finally, the factorization formula © along with our choice © for the jet function lends itself 
to an even more suggestive form for practical applications, namely, as a direct relation between 
the massless and the massive amplitude, 5Vf [P]'( w=0 ) and fM [pM m ) 5 for any given physical process. 
To that end, we exploit the full predictive power of Eq. (@]) and derive the remarkably simple and 
suggestive relation 

M^Uh}^,a^)^ = (7) 

which is the first main result of this paper. 

We have suppressed the color indices in Eq. ©. As we see, the massless amplitude <M [p]'(" ! =°) 
and its massive analogue M W'W in the small mass limit m 2 <C Q 2 are multiplicatively related by 
a universal function zH°). This result is consistent with Ref. [2] (see Section [5J for the detailed 
comparison). The function z( m l°) is process independent and can be viewed as a sort of renormal- 
ization constant (or rather a constant relating two different regularization schemes). This relation 
can be used to predict any massive amplitude from the known massless one, the latter being much 
easier to compute in practice. Moreover, Eq. © includes not only the singular terms in the massive 
amplitude but extends even to the constant contributions (i.e. the mass -independent terms). 

With Eq. © defining the jet function, the function Z^ m ^ is given in terms of the respective 
form factors, 

Z H») ( Ja s ,e) = ,« (J , Ja„,e) (,« (J ,o,«.,e))"' , (8) 

where the index i denotes the (massive) parton and a s is evaluated at the scale /j 2 . Eq. ® explic- 
itly demonstrates the process -independence of the factor Z^°\ While both the massive and the 
massless form factors are functions of the process-dependent scale Q, this dependence cancels in 
their ratio leaving in the factor z( m l°) only the ratio of process-independent scales /u 2 /m 2 . 

Although Eq. ® is valid in a more general setting, and in particular through any perturbative 
order, we will restrict in the following our attention to QCD amplitudes and in particular to those 
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II hi Ih hh 

Figure 1: Feynman diagrams contributing to the vertex y*qq as examples for the flavor classes //, 
hi, Ih and hh discussed in the text. Curly lines denote gluons, double straight lines quarks of mass 
m and single straight lines massless quarks. 



with massive quarks. For this case we will present explicit results for Z}-, up to two loops in 
Section SI Applications of Eq. ([7]) will be presented in Section [5] 

Eqs. © and © are in addition subject to the following clarifications and qualifications. First 
of all, the form factors entering in Eq. ([8]) for Z^"' " 1 are to be understood as being the form factors 
in a theory with either n, + 1 massless quark flavors or n, massless flavors and one heavy quark, 
respectively. In both cases we have the same total number of flavors n* = n { + 1 . Secondly, 
our approach of relating the large logarithms in m to quasi- collinear momentum regions requires 
external massive legs. More precisely, we may define flavor classes, according to the total number 
of heavy quark lines in an amplitude at a given order of perturbation theory. In addition, the flavor 
classes distinguish for a given number of massive lines whether the latter represent external legs or 
form closed internal loops. For the form factor up to two loops, we illustrate the various cases //, 
hi, Ih and hh in Fig.Q] At tree level and one loop, we only have the pair of classes 11 and hi, while 
from two loops onwards we also have the pair Ih and hh. Both these pairs give rise to separate 
relations in Eq. ®. Beyond two loops, yet new flavor classes can appear, see e.g. Ref. [25]. In fact 
a related discussion of this issue has already emerged in the literature during the calculation of the 
NLO QCD corrections to the three jet rate with massive quarks in electron-positron annihilation, 
e + e~ — > qqX, see e.g. Refs. [26,27]. It is also clear how to generalize the definition of flavor 
classes to other types of colored heavy particles such as gluinos. 

Let us finish this Section by pointing out another property of Eq. 0. It is a standard textbook 
knowledge that the two infrared regularizations of any one-loop QCD amplitude, either with a 
quark mass or dimensionally, are related to each other as follows 

ln(m) — > - + finite terms in £ . (9) 

Based on Eqs. © and CD), we conclude in this paper that the proper generalization of Eq. © 
beyond one loop is in the sense of process independent factorization. The factor z( m l°) in Eq. ([8]) 
is invertible and defines the building block of proportionality to all orders in the strong coupling. 
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3 The Sudakov form factor in QCD 



In the previous Section we have presented a factorization that describes the singularity structure 
of QCD amplitudes both in the massless case and in the limit of small masses m 2 <C Q 2 . This 
factorization is valid through any perturbative order and we have emphasized the central role of 
the form factor J W, which specifically includes the QCD corrections. Therefore, in this Section 
we want to focus on J H an d address the issue of its exponentiation. 

To be precise we restrict the discussion here to f w for the vertex j*qq of a photon and an 
external quark-anti-quark pair, i.e. to massive partons in the fundamental representation of the 
SU(/V) -gauge group. Furthermore we confine ourselves to the case of one (heavy) external quark 
line and no internal massive loops, which means we consider the flavor classes // and hi (see 
Fig. [[]). We briefly comment on classes Ih and hh at the end of this Section. The gluon form 
factor f t g l for all-massless partons in the adjoint representation on the other hand, which describes 
the vertex §gg of a scalar Higgs and two gluons is well known, see e.g. Ref. [28]. Also the 
necessary modifications to account for massive partons in the adjoint representation such as gluinos 
in supersymmetric QCD have been worked out to one loop in Ref. [2]. 

Given a photon of virtuality Q 2 (we take space-like q 2 = — Q 2 < throughout this Section) the 
general expression for the vertex function reads 

Tp{k h k 2 ) = ie q u(h) ^y^} q \Q 2 1 m 2 1 a s ) + ^a^q v ^ q \Q\m 2 1 a s )^ u(k 2 ). (10) 

Here the external quark (anti-quark) of momentum k\ (k 2 ) is on- shell with m denoting its mass and 
e q its charge, thus k\ = m 2 (and k\ = m 2 ). The scalar functions and on the right-hand side 
are the space-like quark form factors, which can be calculated order by order in the strong coupling 
constant a s . Results for the perturbative QCD corrections to !f\ in Eq. (flOl) are known through 
three loops in the massless on-shell case [28-31], while the case of on-shell heavy quarks through 
two loops has been considered in series of papers [32-34]. ff" and are gauge invariant, but 
divergent and in dimensional regularization with d = 4 — 2e these divergences show up as poles 
£~^". As we are concerned with the small mass limit m 2 <C Q 2 , we will in the following mainly 
consider the pure vector- like form factor j/ 9 ', since y 2 vanishes for massless quarks. In the 
remainder we drop all indices and define J = . 

The universality of soft and collinear radiation leads on quite general grounds to an expo- 
nentiation of the respective singular terms in the form factor, be it poles in £ or large logarithms 
ln(m) of Sudakov type. This has been well studied in the literature in various approaches [35-39]. 
Moreover, in the massless case explicit formulae have been given up to the next-to-next-to-next- 
to-leading contributions [30,40,41]. However, to the best of our knowledge, an equally valid ex- 
ponentiated representation for the massive form factor in dimensional regularization, which holds 
beyond the leading contributions has still been lacking. In this paper we present it for the first 
time. In doing so we use two complementary derivations based on evolution equations [35] and on 
inclusive partonic cross-sections [39]. 
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Let us start with the former method and recall the evolution equations for the form factor [35] 
2 d , (Q 2 m 2 \ 1 T fm 2 \ 1 /Q 2 \ 

The key input [42] from QCD factorization are the dependence on the hard scale Q which rests 
entirely in the function G and, to logarithmic accuracy, the separation of the mass dependence in 
the function K. Both functions, G and K, are subject to renormalization group equations [35,42], 

l?-Ti G = - limA/ 2 i^fe,a s ,£^) =A(a s ), (12) 



d/j 2 \/j 21 J d/J 2 \ /J 2 

where we assume a s = a s (/v 2 ). Under renormalization group flow both G and K are governed 
by the same anomalous dimension A, because their sum is an invariant of the renormalization 
group. The anomalous dimension A is well known for instance as the coefficient of the 1/(1 — 
contribution to the diagonal splitting functions or alternatively as the anomalous dimension of a 
Wilson line with a cusp [43]. Its power expansion in the strong coupling is currently known up 
to three loops [44, 45] and we use the convention (also employed for all other expansions in a s 
throughout this article) 

A ^) = y L[7t) Ai = A, , (13) 

i=\ 471 i=i 

where we have introduced the shorthand notation a s (/u 2 ) = (X s (/j 2 ) / (4tt) and similarly for the d- 
dimensional coupling to be defined below. For later reference, we also mention that we choose 
the MS-scheme for the coupling constant renormalization. The heavy mass m on the other hand 
is always taken to be the pole mass, thus the renormalization of m imposes the on-shell condition. 
We explicitly relate the bare (unrenormalized) coupling to the renormalized coupling oc s by 

a^5 e = Z as a s , (14) 



where the renormalization constant Z Us in the MS-scheme is given by 

£ ' Ve 2 2 £ J s Ve 3 6 £ 2 ' 3 £ , 

and the bare expansion parameter is normalized as a% = / (4tt) . For simplicity, we always set 
the ubiquitous factor 5 e = (47t) e exp(— e^e) = 1- 

In Eq. (fT2l) all dependence on the infrared sector of the theory, i.e. the structure of the singular- 
ities is described by the function K. The function G, on the other hand, includes all dependence on 
the hard scale Q 2 and is finite for £ — > 0. It is straight forward to solve the evolution equation (fT2l) 
for G. Integration gives 

G(^,a s ,£j = -G(a(<2 2 ,£))- J jA^.e)), (16) 
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with the boundary condition G(a) to be derived by matching to fixed-order results for the form 
factor. 



Working in J-dimensions the solution for G in Eq. (1161) naturally depends on the J-dimensional 
running coupling d(Q 2 ,e). The latter can be expressed as a power series in the usual strong cou- 
pling constant OL s (/j 2 ) evaluated at a scale fj 2 . This relation is now known through NNLO accu- 
racy [30], 

'k 2 \ £ „ 

-2) a(k 2 ,£) = (17) 



+ I lnX Q (3 + Y) ( 1 - X) - 1 - Y - V) } + o {at) , 

which is consistent with the [3-function in J-dimensions [41,46]. Here we have used a s = a s (/j 2 ), 
the obvious boundary condition a(/j 2 ,£) = a s (/j 2 ) and the abbreviations 

$offk 2 Y £ \ v e(l-X) 



x = 1 -«t{[j?) ~ l )> Y = ^r- (18) 

The differential equation for K in Eq. (PT21) is similar to the one for the function G with the obvious 
difference that the scale Q is replaced by the mass m. The solution reads: 

W^-,cc s ,e) = -K(a(m 2 ,e))+ J yA(«(V,e)). (19) 



Combining the above results we obtain the solution of the evolution equation (fTTI) for the form 
factor J , 

fO 2 m 2 \ 

In J ,a„£ = (20) 



■1 / |{G(^ 2 ,E))+/^A(a(V,e))} 







-\ J |{*(a(^ 2 , £ ))-/ ^A(a(V,e))} 



which satisfies the boundary condition f (0,0, 0C s ,e) = 1 and Eq. ([TT]) . Finally, we rearrange the 
above result as, 

(O 2 m 2 \ 

In J ( ^-,^-,a s ,eJ = (21) 

e 2 /^ 2 4 
-\ J ^G(d(^ 2 ,e))+K(d(^ 2 m 2 /Q 2 ,e))+ J ^A(a(fy 2 ,e)) 

W/Q 2 
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Upon expansion of the d-dimensional coupling according to Eq. (flTT) . ln^F develops per power 
of a s double logarithms of Q 2 /m 2 and single poles in £, which are generated by the two integra- 
tions. To be specific, the single poles are governed by the function K in Eq. ([2TI) and are generated 
only by the outer ^-integration. On the other hand, the inner X-integration over A gives only rise to 
logarithms (as long as the infrared cutoff is set by the heavy quark mass). All quantities in Eq. (1211) 
are expressed in terms of the d-dimensional coupling a. In this way all integrations are regulated 
and no singularities other than poles in £ arise. 

Finally, after multiplying with a hard function C, we are in a position to write down the expo- 
nential for the complete massive form factor, which is our second main result, 



with all quantities on the right hand side being functions of a in d-dimensions. Besides the known 
anomalous dimension A [44,45] all other functions G, K and C can be determined in a finite-order 
expansion. This will be accomplished in Sec. @J Before doing so however, there is one more 
feature of Eq. (1221) which deserves some comment. 

As is well known, for renormalization schemes used in QCD based on dimensional regular- 
ization in the MS-scheme, the Appelquist-Carazzone decoupling theorem [23] does not hold true 
in its naive sense. In a theory with n { light and n h heavy flavors (thus «y = rij + n h for the total 
number of flavors) the contributions of a heavy quark of mass m to the Green functions of gluons 
and light quarks expressed in terms of the renormalized parameters of the full theory do not exhibit 
the expected 1/m suppression. The reason here is that the [3-function governing the running of the 
strong coupling constant oc s does not depend on any masses. Neither do the anomalous dimensions 
describing the renormalization scale dependence of all other parameters of the theory. Rather, 
they exhibit discontinuities at the flavor thresholds, which are controlled by so-called decoupling 
constants. 

In the exponential expression Eq. (1221) for the form factor we have used the standard MS cou- 
pling running with n, light flavors. In order to compare Eq. (l22l) or rather its expanded version to 
the fixed-order calculations [32] of the massive form factor, which also employ the MS-scheme, 
but a running coupling with a total number of flavors = n { + 1 , one has to apply the decoupling 
relations. The necessary decoupling constant for a s at flavor thresholds is known to (o^) [47^9] 
(see also Ref. [50]). To relate the two results, that is the expansion of Eq. (1221) on the one hand and 
the perturbative QCD corrections for the form factor through two-loops [32] on the other, we use 
the following relation for a s , 




J ^A(a(V,e))} 



(22) 



K) 




(23) 
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- \ (5C A + 3C F )L,, e + fc A -^C F }( a^) \o({ ^ 



(rij) (if) 

where a s is the standard MS coupling for n { quark flavors expanded in terms of a s for = 
n { + 1 flavors, both evaluated at the scale /j 2 . Eq. (T2~3l uses the pole-mass m. The abbreviation 
denotes 

Eq. (|23l) is correct to NNLO and consistent with the standard |3-function in J-dimensions for all 
terms proportional to L„ e . For the constant terms at al (i.e. those independent of L„ e ) it is accurate 
up to terms of order £. Eq. (T2~3l is to be inserted on the right hand side of Eq. (TTvT) to decouple 
the heavy quark in the J-dimensional coupling. Beyond one loop this generates in particular the 
correct scale dependence to the accuracy required in Section |4] 

Before moving on, we would like to discuss the exponentiation of the massive form factor in 
Eq. (|2TI) from a different perspective. As announced above, our starting point here is the obser- 
vation that in sufficiently inclusive cross-sections, infrared singularities cancel between real and 
virtual diagrams. A suitable example for our purpose is the partonic cross-section of inclusive 
deep-inelastic scattering (DIS) of a massive quark. The purely virtual contributions to this partonic 
observable coincide with the squared massive space-like form factor. 

To extract the form factor, we first derive the all-order exponentiation of the soft singularities of 
the cross-section for the scattering of a massive quark q off a virtual boson V*, i.e. q + V* —^q+X. 
To that end we follow the by-now standard methods for exponentiating inclusive partonic cross- 
sections in Mellin TV-space, see e.g. Refs. [37,39,51-53]. Working in the eikonal approximation 
we obtain 

ln(a(/Y,a s )) = J dxi — jg(l-x,a s ), (25) 


where the function g contains the powers of logarithms ln(l —x) at higher orders of a s and x is a 
kinematical variable related to the Bjorken variable xb, and to be specified below. 

Secondly, we use the fact that the purely virtual diagrams exhibit a simple x-dependence pro- 
portional to 8(1 —x), i.e. in ./V-space they contribute an TV-independent factor. Thus, working in the 
eikonal approximation one can identify the contribution from the squared form factor (to all orders 
in the strong coupling) with the term " — 1" in the factor (x^ -1 — 1) in Eq. (|25l ). The complemen- 
tary, TV-dependent factor is entirely related to real emission diagrams. This way one can identify 
the logarithm of the form factor with the function 

1 fdz 



2 




JjS(z,a s ), (26) 



where z = 1 — x. 
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As it stands Eq. (1261) is not well defined. The reason is that it contains unregulated soft singular- 
ities. Their appearance is not unexpected, since the factor (x 1 ^^ 1 — 1) in Eq. (1251) is constructed such 
that it ensures the cancellation between the soft singularities from the real and virtual corrections. 
Moreover, it is precisely this cancellation that leads to the appearance of the large distributions 
[\n k (l —x)/(l —x)] + (or large logarithms ln(/V) in Mellin space) in the function g. Therefore, 
if one removes the real emission contributions in Eq. (1251) . one can no longer rely on the deli- 
cate balance between real and virtual contributions to regularize the soft singularities. Clearly an 
alternative regularization of the latter is needed to render Eq. (1261) meaningful. 

Since in this paper we are interested in regularizing the soft divergences in the massive form 
factor (or in any other amplitude) dimensionally, and in line with our previous discussion, we 
modify Eq. (l26l) by replacing the usual coupling a s with the J-dimensional one a as defined in 
Eq. CEZ]), 

l 

ln(5(a,)) = ~/yg(z,a). (27) 
o 

We stress that the function g in Eq. (1271) is the same one that appears in the cross-section in Eq. (l25l) . 
The effect of the d-dimensional coupling is rather transparent, as it supplies additional powers of 
the factor z~ e , see e.g. the left hand side of Eq. (fTTT) . which allows to regulate the z-integration in 
Eq. d27j in the limit z -> 0. 

For the derivation of the required hard cross-section for the process q + V* q+X we directly 
build on previous work on the exponentiation of massive cross-sections at next-to-leading logarith- 
mic accuracy [54,55], where the light quark initiated process qi + V* — * q+X was studied. Since in 
this work we are interested in the corresponding process initiated by heavy quarks, q + V* — > q +X, 
one has to modify the analysis of Ref. [54]. One possible option is to repeat the considerations of 
that reference keeping a non-vanishing mass for the initial state quark. However, a much simpler 
alternative is to express the coefficient function for the ^-initiated process as a convolution of a 
perturbative distribution function for the initial-state heavy quark q and the coefficient function for 
the process q\ + V* — > q+X both evaluated at a common factorization scale /jf. Since we are inter- 
ested only in contributions that are enhanced in the soft limit and suppress power corrections with 
the quark mass m, only the q — > q component of this distribution function is required. Moreover in 
the soft limit this function with space-like kinematics coincides with its time-like counterpart (see 
e.g. Ref. [56]). All components of the time-like perturbative fragmentation function D are known 
through two loops and can be found in Refs. [57, 58]. 

The exponential structure in the soft limit of the perturbative fragmentation function D of a 
heavy quark [19] is well understood [56,59,60]. In Mellin-/V space we have 

I i^-i-i f *f dk 2 1 

ln(D(N)) = J dx | H (cc s ((l -x) 2 m 2 )) + J -^A (a s (k 2 )) J , (28) 

(\-x)hn 2 

where the anomalous dimension A is the same as the one appearing in Eq. (fT3l) while H is a new 
function. 
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The exponentiation of the coefficient function for the process qi + V* — ► q + X was clarified in 
Ref. [54]. With the same anomalous dimension A and a new function S the result reads, 



2 ijr2 



1 (l-xYM- 

r r^ -1 - I f r dk 2 } 

ln(o M (iV)) = y ^ ^T^T { 5 ^ s((1 _x)2m2) ) + y lir A W^)) } • W 

where in the limit m 2 <C Q 2 the scale M equals M 2 = Q 4 /m 2 . 

Let us briefly recall a few basic facts [54] about the derivation of Eq. (|29l) . The variable x, 
0<x< 1, is therescaled Bjorken variablex= (1+m /Q )xb. The upper limit of the /^-integration 
follows from kinematics and in the center-of-mass frame it is determined from the light quark 
energy E : k 2 <AE 2 {\ — x) 2 . Moreover one has 

Q 2 + m 2 

2E ~ f=== (30) 

v /(l_ J ) 2 2 +m 2 

Since we are working in the soft limit (1 — x) — > 0, it is obvious that in the massive case Eq. (1301) 
leads exactly to the scale M in Eq. <[29l) , while for m = it reduces to the well known expression 
of the massless case [39]. 



Convoluting Eqs. (1281) and (1291) we obtain the desired coefficient function for the sub-process 
q + V* —> q + X in the soft limit. One can see that the dependence on the factorization scale drops 
out as it should. Following the procedure outlined around Eqs. (|26|) . (1271) above, we finally obtain 
the Sudakov exponent for the massive form factor: 

1 z 2 M 2 2 

A F = - 1 -fj! [ S(d(z 2 M\e))+H(a(z 2 m\e))+ J ^- A (a(* 2 ,e)) } . (31) 

z 2 m 2 

The complete form factor is obtained by multiplying the above exponent with a hard function Hp: 

^l/?v ,as,£ ) = Hf<kQ ' 5 ^' E )' e ) ex P{M- (32) 

Here all functions A, H, S and Hf have perturbative expansions analogous to Eq. (fT3l) . They 
can be obtained to a given order in a s by matching for instance to the full calculation for J . In 
addition, independent information on H and S arises also with the help of Eqs. (|28l) or (|29l) from 
the calculation of the perturbative fragmentation function D or the hard partonic light-to-heavy 
DIS cross-section G qi -+ q - 

The hard function Hf has an expansion in £ but is finite in the limit £ — » 0, since all soft 
poles are collected in the exponent A^. To completely define the hard function Hf one has again 
to specify the definition of the coupling a s appearing in Eqs. (I3TT) and (1321) . This is the usual 
MS coupling defined in Eq. (fT4"l) but running only with the number of light flavors n,. The same 
number of flavors appears also in the anomalous dimensions, see e.g. Refs. [56,57]. Thus, in 
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order to compare Eq. (1321) to the fixed-order calculation available in, say Ref. [32], with a coupling 
constant oc^ for = n t + 1 flavors, we again have to apply the decoupling relations [23,47-49] in 
the form of Eq. (l23l . 

A comparison to the exact two-loop calculation of the vector form factor [32] shows that 
Eq. (|32l correctly predicts all soft terms ~ 1/e , k > 1 including their logarithmic mass de- 
pendence, while it does not control the powers of ~ ln^(m) at order £° which are of collinear 
origin. From the viewpoint of the exponentiation of soft singularities these latter logarithms must 
be included in the hard function Hp. 

However, at the same time one expects that all pure collinear logarithms exponentiate as well. 
This feature is unrelated to the soft-gluon exponentiation discussed above but rather to the standard 
parton evolution equations (DGLAP). Here we recall the analysis of Ref. [39] where the effect of 
collinear radiation in the outgoing jet results in modifications of the naive eikonal exponentiation. 
The additional collinear contributions in the final- state are taken into account by constructing a 
DGLAP- like evolution equation for the corresponding jet function. The latter, in turn, contributes 
to the well known DIS anomalous dimension B [37, 39]. In the massive case the virtuality of the 
final state is of order (1 — x)Q 2 + m 2 and does not vanish in the soft limit which brings additional 
ln(m) terms. In this paper we will not elaborate on that point further, as all purely collinear logarith- 
mic terms can be read off from the exponentiated expressions of the form factor given in Eq. (|22l . 
This picture is consistent with fixed-order calculations of the perturbative fragmentation functions. 
Indeed, one can easily verify that the logarithmic contributions in the one-loop form factor f \ for 
E° and £ coincide with the pure virtual contributions to the one-loop fragmentation function D\, 
see for instance Eq. (45) of Ref. [57]. Unfortunately, the two-loop virtual contributions cannot be 
extracted from Ref. [57]. In the next Section we will elaborate on this relation. 

Before completing this Section on the massive form factor, we would like to address the ques- 
tion of its massless limit. Reconsidering Eq. (l20l) . i.e. the exponent of Eq. (|22|) . it is obvious that 
the limit m — > is smooth. The contribution from the infrared function K vanishes for m — > 0, while 
the contribution from the G function is similar to the one in the well known massless case [40,41]. 
The counterterm K in Refs. [40,41] can be expressed as 

l 

/dx 
— A(a(xf?,e)) , (33) 

o 

which is consistent with Eqs. (|20l ) and (1221) after changing the integration boundaries according to 
Refs. [40,41]. Moreover, the function G in the massive case should be related to the G function of 
the massless case up to possibly constant difference. As we will demonstrate to two loops in the 
next Section, the two functions in fact do coincide to all known orders in £, i.e. the massless limit 
of the massive result Eq. (1221) requires setting both m = and C = 1 . 

1 We would like to thank S. Catani for an interesting discussion on this point. 
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Therefore, the resummed quark form factor reads in the massless case 




(34) 




o o 



with the boundary condition J (0,0,a s ,e) = 1. Now In f develops double poles in £ per power 
of a s from the X- and the ^-integration over the anomalous dimension A. In Eq. (l34l) we have 
identified the initial condition G (a) of Eq. (fT6l) with the sum of two functions B + h. The physical 
interpretation of the new functions B and h, which also have expansions in the J-dimensional 
coupling, follows nicely from the previous considerations of inclusive DIS scattering. Following 
Ref. [39] one can identify the function B with the coefficient governing the evolution of those large 
logarithms ln(iV) in inclusive DIS scattering associated with the final state jet function. This has 
recently also been pointed out in Ref. [61]. The function B is known to three-loop accuracy [53] 
from explicit DIS calculations [25,44,45,62]. The new contribution h on the other hand can be 
thought of as the massless limit of the function H in Eq. ([3TI) . To determine h in a perturbative 
expansion we match the above exponent to the known three-loop result for the massless form 
factor [30] and we have checked that the oc^ prediction based on Eq. (1341) agrees with previous 
results in the literature [30]. It is interesting to note that unlike the standard expression for the 
massless form factor as given e.g. in Refs. [30,40,41], the result we propose in this paper is 
comprised of well defined integrals. Moreover, we can directly interpret the respective parts as 
the three-loop contributions of the form factor to the DIS coefficient functions and the splitting 
functions respectively. 

Finally let us briefly comment again on the various flavor classes, since the previous discussion 
was entirely focused on the flavor classes // and hi. Beyond one loop we have for instance the 
contributions to f from the class hh. These contributions are finite after performing the ultraviolet 
renormalization, but they still do contain Sudakov logarithms of the type \n k (Q 2 /m 2 ). In fact, up to 
two loops all remaining large logarithms \n k (Q 2 /m 2 ) in the heavy quark form factor not accounted 
for by Eq. (1221) are entirely related to the self-energy contributions of a heavy quark, i.e. the 
diagram denoted with hh in Fig.Q] It is well known [63,64] that these contributions obey Sudakov 
exponentiation similar to Eq. (1221) . although with different integration boundaries and evaluated 
at the matching scale /j 2 = m 2 . Thus, we can introduce the function J hh which exponentiates the 
logarithms in the flavor class hh, 



Eq. (l35l) is to be evaluated at the scale i? = m 2 and to be restricted to the purely fermionic contri- 
butions with the heavy quark pair coupling to the external boson and one additional virtual heavy 



lnThh(Q 2 ,m 2 ,a(^,e),E) = 



(35) 
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quark line. Eq. (1331) contains the same functions A and K as Eq. (1221) . but a different function G'. 
As a matter of fact, its structure follows directly from integrating Eqs. (ITTT) and (fT6l) under the 
condition that the infrared region is cut off at the scale m 2 . See e.g. Refs. [63,64] for details on the 
finite-order expansion of exponentials like Eq. (1351) . We address this issue in future work. 

4 Fixed-order expansions and resummation coefficients 

In this Section we will give the finite-order expansions for the various quantities, in particular for 
the factor z( ra l°) of Eq. ([8]) and the form factor J of Eq. (|22)) . All formulae in this Section use the 
MS-scheme for the coupling constant, while the heavy mass is always taken to be the pole mass 
(on-shell scheme). In addition, as in the previous Sections, we limit ourselves to the contributions 
in the flavor classes // and hi. Throughout this Section ny denotes the number of massless flavors. 

4.1 The factor/ '" °) 

Let us start with the perturbative QCD expansion of Eq. ([8]) for the factor z( m l°) which we defined 
as the ratio of the massive and the massless form factors for a given parton and which we write as 
an expansion in terms of the renormalized coupling a s (/u 2 ) = CL s (/j 2 ) / (4k): 

^ 0) (^o,,e) =l + EW i z[j ) . (36) 

The quark form factor in massless QCD is known to three loops [28, 30], while the form factor 
of a massive quark is known for arbitrary values of the quark mass through two loops [32]. The 
expansion coefficients (|36"1) in the case of a (heavy) quark q read 

4,1 = C f {A + ^ +L J +Lf , + 4 + C 2+E (| + | + (4 + C 2 )L f , + 8 + |-|c3) (37) 

- 2 (i4 + ( 2 4K + ( 8 4-^K +i6+ ^-! + ^ 2 )} 

+ 0(£ 3 ), 

Zg> = C/^ + ^{c F 2 (4L, + 2)- 1 -±C F C A +n f C F ] (38) 
+ ^ {Cf 2 + *Lp + y + 2&) + C F C A (- y L, + y - fc) + nf F (Jl„ - } 

+ i K + < + ( 17 + 4 ^ )L " + T " 4 ^ 2 + T^ 3 ) 

^ 2/ 4 r4 8 r3 tin , r N r 2 / 83 64 „ \ 561 61„ 22 „ 

+ C F 2 (y^ + -Lj + (17 + 4^)L 2 + (y - 8^2 + y CaJ^ + -y + y C2 - y ^3 
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_ 48 l„2fe - lltf) + C F C A {^Ll + (f - 2fc) L> + (If + f - 30?,) Z, 
12877 323, 89, ,„ „, 47, 2 

+ w + Tr fc+ ^ 3+241n2fc -y fc 

„ C 2,3 13,2 ( 17 8,\, 1541 37, 26, \ 
+ "/CH"9 L J " T L 2 + ("27 " 3 & K " 324~ T & " + ° (e) ' 

where two-loop contributions arising from virtual heavy flavor lines are omitted (see Fig. Q]) and 

Z, = l„g). (39) 

In presence of heavy flavors the form factor of a massless quark gets mass-dependent contribu- 
tions at two loops from the diagram Ih in Fig. [TJ We will not consider such two-loop contributions 
in this paper. Unlike the massless quark, however, the gluon form factor receives mass -dependent 
corrections starting from one loop. These have their origin in the one-loop heavy flavor insertion in 
the tree-level gluon form factor. It is clear that for the gluon form factor the classification of Fig. Q] 
has to be suitably adapted by counting the number of (internal) heavy lines. It is easy to work out 
the one-loop result for the gluon z( m l°) -factor in Eq. (l36l) (see the Appendix for the all-orders in 8 
result) and it reads: 

4t - -{-|-^ + <-^-T) +E2 (-^-|^ + f )}+°< e3 >- (40) 

In addition, the following comments on Eqs. (1371) , (1381 ) above are in order. First of all, we 
have to supply the £-expansion of the massive form factor including terms of order £ 2 at one loop, 
because of the singularity structure with l/£ 2 -poles in massless one-loop form factor. Since the 
(£ 2 ) term of the one-loop massive form factor is not available in the literature we have calculated 
it following the setup of Ref. [32]. Details are given in Appendix lAl One can easily verify that this 
term produces a finite contribution at two loops for an amplitude with n h external massive quarks, 
fil external massless quarks and n g external gluons, which would be proportional to nflp 1 + n g C F C A 
times the Born term. 

Secondly, there is one important detail about the scheme for definition of coupling constant and 
masses. We assume the pole-mass definition for the heavy quark mass m as well as the standard 
MS coupling defined in Eq. (fl4l) . Note that this definition for the coupling differs from the one used 
in e.g. in Ref. [32] (and other references on higher order corrections for massive processes) where 
the coupling renormalization includes also the factor T(l +£)exp(£^). For consistency with the 
massless calculations, we have performed a finite renormalization of the result in Ref. [32]. The 
necessary relation is given by 

= aJl+a s -(p -l) ( r( , 1 + e) -1^+0(^)1, (41) 
Ref. [32] s \ S £V H 3 J Vexp(-£YE) J V s) ) ' 

where we put the factor (47t) e exp(— £Je) = 1 for simplicity and Po = 1 — 2/3fly. It is easy to 
see that through two loops this amounts to the following finite correction (see Eq. (|43T) below for 
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definitions of 7i) to the results presented in Ref. [32] 



7i 



_ = 7i 

MS 



Ref. [32] 



+ ^f(fe 2 + o(e3)) 5 , 



(42) 

Ref. [32] 



(ra|0) 

Finally we would like to elaborate on the relation between the factor Z and the heavy 
quark perturbative fragmentation function we discussed in Section [3] preceding Eq. (l34l) . At one 
loop, the virtual contribution to the fragmentation function was explicitly calculated in Ref. [57] 
to all orders in £. We present this result in Appendix [A] in a particular form prior to collinear 
factorization and one can easily verify by a direct comparison that its expansion through (e 2 ) 
coincides with the factor zPj . Moreover, Ref. [57] also contains the purely virtual fermionic 
contributions (i.e. proportional to the number of light flavors) at two loops. In terms of the usual 
renormalized coupling we have found the former to be in agreement with the terms proportional to 
rij: in the function Z>J to all powers in £ appearing in Eq. (1381) . This observation indicates that the 

factor zf"^ of Eq. ([8]) indeed coincides with the virtual corrections to the collinearly unfactorized 

perturbative fragmentation function and one may actually view the complete agreement between 

(m\0) 

all known terms of the two functions as a check on the derivation of Z r , 1 . 

m 

Although the latter object is not known to the level we have presented here for the function 
Z^ 1 ' the apparent coincidence allows for an interesting alternative interpretation of that function 
by relating it to the field renormalization constant of a heavy quark in light cone gauge n ■ A = 0. 
Indeed, in the approach of Ref. [57] to calculate the fragmentation function, the purely virtual 
corrections are nothing but insertions of self-energy type in external on-shell legs in this particular 
gauge. Clearly, it will be very interesting to further develop this line of reasoning. 



4.2 The form factor J 

Next, we want to perform the finite-order expansion and matching of Eq. (1221) for the heavy quark 
form factor <J . Subsequently, with all functions G, K and F determined we will then be using 
Eq. (1221) for predictions of perturbative results at higher orders and derive explicit results at three 
loops. To that end we perform the integrations in Eq. (1211) after inserting the perturbative expan- 
sions of all quantities and simply evaluated resulting integrals. Details on this procedure may be 
found in Refs. [30,41]. 

For the (ultraviolet) renormalized massive form factor with space-like virtuality q 2 = — Q 2 < 
and in terms of the renormalized coupling <X s (fj 2 ) we have, 

7 (^,^Mf)A = 1 + EW>, (43) 

With the convention of Eq. (TT3l for the expansion of A, G, K and C and setting the scale to /u 2 — m 2 , 
we find 

1 rl 1 . ,i 1 A ~ 1 r 1 A , 1 ,i 

7\ = -{-A l L+-(G l +K l )j--A l L 2 --G l L + C l +e[—A 1 L 3 + -G 1 L 2 j (44) 
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- £2 {^ AiL4+ i^ GiL3 } + 0(£3) ' 

^ { l -A\ L 2 + l -A i (Gi + *i - po)^ + i (Gi + *i ) (Gi + K X - 2f3 ) } (45) 
+ ^-^A 2 L 3 -^A l (3Gi+K l )L 2 + ^(A 2 -G 2 -K l G l +2A l Ci)L+^(G2 + K 2 ) 
+ -C X (Gi + Ki ) } + ^rA 2 L 4 + ^A i (7Gi + K x + 2(3 )L 3 + ^Gi (2Gi + #i + 2(3 )L 2 

- ^(A 2 +AiCi)L 2 - ^(Ga + GiCOL + Cs + ej-^AfL 5 - ^A X {\5G X +K X +6(3 )L 4 

- ^Gi(4Gi + 6p )i 3 + ^(2A 2 +A!Ci)L 3 + ^(2G 2 + G l C l )L 2 } + o(e 2 ) , 

i{-U?L 3 + ^A 2 (Gi + ^ -2p )L 2 + 4 A i( G i + *i)(Gi +*i -4Po)^ (46) 
+ l -A$ 2 L+ i(Gi +ZO(Gi + *i -2Po)(Gi + *i -4p )} + ^{-1a 3 L 4 

- ^A 2 (2G! + *i - (3 )L 3 + ^Ai(A 2 +A l C l )L 2 + -UiPo(3Gi +K l )L 2 
Id o Id 

-^A 1 (G 1 +K 1 )(5G 1 +K 1 )L 2 + ^A 2 (3G l + 3K 1 -4^ )L+^A 1 (3G 2 + 3K 2 -4^ 1 )L 

- ttGi(Gi +£i)(Gi +£i - 2p )L+ 7 A!d(Gi +Zi - p )L+ ^Ci(Gi + £i) 2 
Id 4 o 

+ A. (Gi + (3G 2 + 3^ 2 - 6(3 C! - 4Pj ) - ^ (G 2 + * 2 ) } + i { ^A 3 L 5 

+ -i-A 2 (25G 1 + IKi + 4(3 )L 4 + ±-A x ( 19G 2 + K 2 + UK X G X )L 3 
iyz yt> 

+ T? A iPo(4Gi 4-^)L 3 - -Lai(3A 2 + 2A 1 Ci)L 3 + i-Gi(Gi +#i)(3Gi +^i + 2|3o)L 2 
4o Id J2 

- ^A 2 (2Gi +Ki)L 2 - ^ 2 A l (5G 2 +K 2 + 6G1C1 +2^Ci)L 2 + \-A x ( K 32C A - \35C F )L 
8 16 3d 

+ ^-(2A 3 + 3A 2 C 1 +6A 1 C 2 )L+l(-3G 1 G 2 -2K 1 G 2 -K 2 G 1 )L-]-G 1 C 1 (G l +K 1 )L 

12 o 4 

+ ^(G 3 +^ 3 ) + ^C 1 (G 2 + ^ 2 ) + ^(G 1 +^ 1 )(18C 2 + 32C A -135C F )} 

- ^-A]L 6 - ^tA 2 (6Gi + 3|3o)L 5 + -^rAi(16A 2 + 7AiCi)L 4 - J-A^G 2 
04 o4 yo Jo4 

+ 30AiGi + £ 2 + 90p Gi + lOpo^i + 16p^)L 4 + — A 2 (13Gi +4i^i +8(3 )L 3 

48 

+ ^-A 1 (19G 2 +^ 2 + 4(3 1 )L 3 + i-A 1 C 1 (7G 1 +^ 1 +2p )^ 3 -^G 1 (9G 2 +^ 2 
4o 24 yo 

+ %K X G X +22p Gi + lOpo^i + 16po)^ 3 " 7(^3 +A 2 Ci)L 2 + -^Gi (^ + 4(30^ 

4 16 

- -U 1 (18C 2 + 32C 4 - 135C F )L 2 + — G 2 (9Gi +4Ki + S$ )L 2 
72 Id 
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+ IgiCi (2Gi + + 2p )L 2 - l - (G 3 + G 2 d + GiC 2 )L + C 3 + O (e) , 
where again contributions arising from virtual heavy flavor lines are omitted (class hh) and 

*■= >■>(!)■ <«* 

In the quantities G3 and C3 in Eq. (1461) we have also absorbed all constant contributions from the 
decoupling relation (1231) at order a^. All these terms are independent of and can potentially 
include contributions of order £ at aj which we did not write out explicitly in Eq. d23l . Results for 
fx at a general scale fj 2 7^ m 2 can be derived from Eqs. (I44~l)-(l46l) by standard methods^]. They will 
be presented elsewhere. 



Explicit results for Ji in Eqs. (144)) -(|461) can be obtained with the help of the known coefficients 
of the cusp anomalous dimension A(a s ) due to Refs. [44,45,65], 

Ai = 4C F , (48) 
268 „„ \ „ / 40 \ 



A2 = C F C A (— -S^ +rif C F (--), (49) 

A3 = C F C/(f -^ 2 + f, 3 + fc 2 2 ) (50) 
+»/C,C A (-f + ^ C2 _^ 3 ) + n / C/(-^ + 32C3) + n/C F (-^). 



The respective coefficients for G(a s ) and ^T(fl s ) read, 

/ 28 \ 

Gi = -6C F +£C F (-16 + 2C 2 )+£ 2 C F (^-32 + 3^ 2 + y^J (51) 
+ £ 3 C F (-64 + 8^ 2 + 14^3 + ^2 2 ) + (£ 4 ) , (52) 
G 2 = C F 2 (-3 + 24^ 2 -48^)+C F Q(-^-^ 2 + 52^ 3 ) (53) 

*l = 2C F , (54) 
K 2 = C F 2 (3-24^ 2 + 48C 3 )+C F C A (— + 30C 2 -60^ 3 )+n / C F (---4^ 2 ). (55) 

Here we have included higher orders of £ in the anomalous dimensions, to ensure that all large log- 
arithms in m are actually generated entirely by the integrations over £, and X in Eq. (I22l) . Although 
this is not a compelling choice it captures all structures in the exponential, which are universally 
related to parton dynamics. This is contrary to "minimal" versions proposed e.g. in Ref. [24]. 



We take the opportunity to point out a typographical mistake in Eq. (62) of Ref. [32]. The following term 
C A C P ±(l - ^ - (7-b)) "(0,-*) should actually read C A C F ^ (§ - ^ - J^)h(0, X ). 
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For the coefficients of the matching function C(a s ) we find, 

d = c F (4 + ^ 2 )+eC F (8 + ^ 2 -^3)+e 2 C F (l6 + 2^-^3 + ^2 2 )+o(e 3 ), (56) 

C 2 = C F 2 (30 + 55k-36C 3 -48^ 2 ln2-^ 2 2 ) (57) 

„ / 2387 71 479 3 „ 2 \ 

+ QC F (-_ + -ta + — Cs +24^ 2 ln2 - -tf) 

„ /356 37,. 38„ \ 

Putting everything together, including the terms of order £ 2 at one loop (see Appendix© we arrive 
at the following results, 

7 X = C F {^(2L-2)-L 2 + 3L-4 + 2C 2 + £^L 3 -^L 2 + (8-^ 2 )L-8 + 2^ 2 + 4^) (58) 
+ e2 (_^ + l t 3_ (4 _l fe)Ll+(l6 _ 3 52 _M b)t 

_ 16 + 6C 2 + yC3 + yC2 2 )} + 0(e 3 ), 

C F 2 1 ^ (2L 2 - 4L + 2) + - (-2L 3 + 8L 2 - ( 14 - 4^ 2 )L + 8 - 4^ 2 ) (59) 
+ 7 -L 4 - ^L 3 + (f -4^ 2 )l 2 - - 32C 3 )L + 46 + 39^ 2 -44^3 -48^ 2 ln2 
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r 1 /2 2\ 1/10 10\ 2 , 19 o /209 4„ \ 212 

+ CpUfl { —L ) H — ( L H ) L H L — ( + -C 2 U + 

F/ l£ 2 V3 3/ e\ 9 9 / 9 9 V 27 3 s / 

_« fc _^ + e( ^_« 1 3 + ( ^ + fc)l2)}+0(le) + 0(E2) , 
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iF 3 = C F 3 {^(^L 3 -4L 2 + 4L-l)+i(-2L 4 + 10L 3 -(22-4C 2 )L 2 (60) 

1/5. « 34, /137 



+ (22 - 8£ 2 )L - 8 + 4^ 2 ) + - (|l 5 - *-L A + (Hi - 6^ 2 )l 3 - (89 - 56^)l 2 

+ (129 + 88^2 - 136^ - 96^ 2 ln2 - - L 6 + yL 5 - - ^ 2 )l' 

494 , 17„ 268 „ \ t3 i , „ r „ (If 22 t7 , 44 r 22 \ , 1 /ll 



+ g-«,)tf-g-|c. + *)L-? + ^ + < 1 ) + I(^ 
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" " ^i)L 3 + (-^- + y C2 - 5Q^)L 2 - - -£ 2 - — & - 48^2 ln2 

4- 26 C 2 V > \ n ^,/4289 16 N 4 / 6260 97 232 \a\ 
^ .rl/242 r 242\ 1 / /2086 44.. \ 1690 44^ 44.. 



+- 

£ 



1 //245 536.. 44.. 88„ 2 \ \ 121 4 /2869 44,. \ 3 i 

„2 r 1 / 4 ,2 8, 4\ 1/ 2 , 4,, /10 4„\, 8 4, 

+ I (-V + « £ , _ (481 10 3 2 ,599 _ 2 8 ^ 1 5 _ S 
£ V 9 9 V 27 3 / V 27 3 s 3 s V / 2 5 
/1016 29„\ r3 i f 1/ 88 r 88\ 1 //668 8„ \_ 

596 8„ 8„ \ 1 / /418 80 „ 56„\\ 22 r4 /974 8„ , 



81 9* 9' / e\ V 81 27" 9' / / 27 V 81 9 
8 \ 1 { 40. 40\ 1 / 8 _\ 2 . 4 76_ 3 



-) — » I L+ — + L L H L 

277 £ 2 V 81 81/ £V 81 / 27 81 



+ c/{ l -[-l5L)}+C P C A { l -{^L)} + O(L°e.- , ) + O(L 2 e ) + o(t). 



Further improvements on the accuracy of the three-loop prediction require an extension of the 
two-loop result J2 to order £. We will return to this issue in a future publication. 

Let us close this Section with a few comments. First of all, it is clear we can obtain a three-loop 
prediction, i.e. the coefficient zj^ in Eq. (1361 ) for the factor Z^"' -* from the exponentiated massive 
form factor in Eq. (1221) with the help of Eq. ( |60l ) and the known three-loop results in the massless 
case [28,30]. Next, putting the discussion in a broader perspective, we note that exponentiations 
similar to Eq. (1221) have also been studied for electroweak interactions [63, 66] in massive gauge 
theories, where large logarithms in the mass of the gauge boson appear. There, the resummation 
has been used as a generating functional for Sudakov logarithms at higher orders. Of course, 
those details of the exponentials which depend on the infrared sector of the theory are modified in 
comparison to Eq. (1221) . However, it is rather striking to observe that the coefficients for G\ and 
G2 from our determination in Eqs. ([FT]) and ( |53l agree precisely with the values for J^ 1 ) , t^ 2 > in 
Ref. [66] (up to an overall factor 1/2 due to different normalizations). In both cases, the relevant 
coefficients control the single logarithm L at the respective order. 



On top of this, it is even more striking, that the very same coefficients G\ and G2 from Eqs. (1511) . 
(l53l for the form factor of a massive quark also coincide with the corresponding results in massless 
case [28,30] (up to an overall sign from different definitions). This observation, which calls for an 
explanation, suggests a universality of the function G which extends even to higher orders in £, see 
e.g. Eq. (I5TT) . It also offers the chance for a conjecture about the coefficient of the single logarithm 
L for 53 in Eq. ( |46l ) purely on the basis of the corresponding massless result, provided, of course, 
that all necessary terms to higher order in £ up to two loops are known. 

The potential consequences of such a universal nature of the function G would be rather in- 
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teresting. For instance, in the massless case there exist additional relations between the functions 
T and J w\ i.e. the form factors for the vertices y*qq and §gg. These relations manifest them- 
selves in underlying structures for the respective function (z = q,g) [28,61,67] such that one 
can decompose the resummation coefficients G^ in the massless case according to 





= 2 






= 2 




of 


= 2 


(4 


where i 


= q,g 


and 



(61) 
(62) 



F = £ [F-F(e = 0)] . (64) 
Here (and only here), the functions B$ (not to be confused with the ones given in Eqs. (|69|)- 



(1711) ) denote the coefficients of term with 8(1 — x) in the n-loop diagonal MS splitting functions 
l ' [44,45], while the universal functions /« exhibit the same maximally non-Abelian color 
structure as the A$ [43] up to the factor C A /C F , i.e., ff = C A /C F f) q \ see Ref. [28] for details. 
With the help of Eqs. (|6T|)-(|631 one could translate the exponentiated form factor in Eq. (I2TI) for 
heavy quarks immediately e.g. to the case of gluinos g with mass mg, all resummation coefficients 
up to three loops being known. 

For completeness, let us finally mention also the functions H and S of Eq. (I3TT) as well as B and 
h from Eq. (l34l) . The function H is already known through two loops from Refs. [56,59], while 
the function S was evaluated in Refs. [54, 55] to one loop. Using Eq. (l32l) and matching it to the 
fixed-order calculation for if 2 we can extract in particular the two-loop coefficient S2 from the term 
a^/eL°. The explicit results for H(a s ) of the massive form factor in Eq. (1321 read, 



Hi = -4C F , (65) 
.(^ + 8C2 -36C 3 ) + n / c4 



H 2 = C F C A ( — + SC,2-36t >3 )+n f C F —, (66) 



and for S(a s ) 



51 = -4C F) (67) 

52 = C F C A (-^ + ^ 2 + 20b)+n f C F ^. (68) 

In order to have a self-contained presentation, we also give the perturbative expansions of the 
coefficients B(a s ), 

Bi = -3C F , (69) 

0/ 3 \ / 3155 44 n 

Bi = C F 2 (~ + 12fc- 24^3) +C F C A (-^ + -^ +40^ 

+ C ^/(^-^ 2 )' (70) 
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B 3 = C F 3 (-y-18C 2 -68C 3 -^2 2 + 32C 2 C3 + 240C 5 ) 

,/ 712 272 9 \ 

+ C A C F 2 (-46 + 287^2 - — Cs - — ^2 - 16^2^3 - 120C 5 ) 



„ ~ ( 599375 32126. 21032. 652. 2 176. . 
+ Q 2 C F (-_ + _ ^2 + & - ^ ~ ~ 232^ 5 

n 2 /5501 cn . 32 r \ ^ 2 / 8714 232,, 32 „ \ 
+ C A ("sT " ^ + 9"W +C W ("129 + 27^ " 27 W 
/ 160906 9920 776 208 2 \ 

and for /i(a s ) 

fci = -3C F + £C F (-16 + 2C2)+e 2 C F (-32 + 3^2 + y^)+e 3 C F (-64 + 8^ 2 (72) 
+ 14C3 + ^C F (-128 + 16^2 + f Cs + ^tf - f + f Cs) 

+ e 5 C F (-256 + 32^2 + f C 3 + £tf - 7^3 + f Cs + ^ 3 - ?C3 2 ) , 

&2 = C F 2 (-^ + 12^2-24^) +C F C A (-^-yfc + 12? 3 )+ii / C f (y + yfc) (73) 

^2^ 1 176 „ 2 \ ^ ^ / 70165 575,, 520. 
+ £C F 2 (-- + II6C2 - 120^3 - — C2 2 ) + eC F C A — £ 2 + — 

176. 2 \ „ /5813 74. 16. \ . 2 /109 .„„, 432. 2 

+ -j-tf) + enj€ P + — C2 - yCs j + e 2 C F 2 (— + 437£ 2 - 736£ 3 - — C2 2 

y AO „\ 2^ ^ ( 1547797 7297. 24958^ 653 , 2 356 „ „ 
+ 1 12C2C3 - 48^ 5 ) + e 2 C F C A ( — — C2 + ^7-^3 + -^-C 2 2 - — C2C3 

.... \ 2 „ / 129389 850,, 1204^ 7. 2 \ , 2 /1287 2991. 
+ 204^ 5 ) Wn f C F (_ + — fc - — ^3 - 3 C2 2 ) +B 3 C F 2 (_ + — C 2 

- 3614^ - 508^ + 104^2^3 - 72^ 5 + + 1072^ 3 2 ) + e 3 C F C A {~^- 

155701. 308810. 100907. 2 478. . OAn „ 1618. 3 2276 . 2 

i«- &H — 81— & + ^l80- 5a --^3 + 840C 5 -^ 2 -^3 

3 ^ / 2628821 8405. 16340. 1873^ 2 44, , , . \ 

h 3 = C F 3 (-y-18C 2 -68C 3 -^C2 2 + 32C2C3 + 240^)+C F 2 C A (-^ (74) 

1235 2296 856 2 lfir r „ nr \, rr 2 / 16540 22286 
+ C2 3— C, 3 + -jyC2 - 16^2^3 - 120^5 J +C F C A ^ — ^7~C2 



1344. 2 \ 2 /^y 14t)„ 4UU 2 

7 F A V 27 27 b 3 b 15 b / 7 F V 27 27 b 



where the function B is known to three loops from Refs. [37, 39, 53, 62], while the function h has 
been derived by matching Eq. (|34"1) to the respective fixed-order calculation starting from the single 
pole terms a"/e. 
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5 Applications 



Here we want to demonstrate how the previous considerations can be applied to derive the structure 
of the singularities and all large Sudakov logarithms in higher order QCD corrections to partonic 
scattering processes. Let us start with the general 2 — > n scattering processes of partons p, in Eq. (QQ) 
and consider Eq. © in a perturbative expansion in oc s . We want to present the explicit relations 
between corresponding amplitudes with and without parton masses {m,}, in our notation \ M p s m .x) 
and |5Wp,{ m ,.=o})- Throughout this Section, we consider (ultraviolet) renormalized quantities and 
we define 

CO 

I^p) = IW>p ( °), (75) 

i'=0 

and any overall powers of a s typical say, for jet cross-sections at hadron colliders, have been 
absorbed in the notation. We can then express Eq. © for a general process (OQ) in an expansion to 
second order in a s as, 

\M {Q ) ,) = nx ), (76) 

I p,{m,}/ I p,{m,=0}/ ' v ' 



) - - V Z (1 W (0) ) + ) (77) 



e {all legs} 



;) = 1 E (zff- 1 ^) 2 )^ m ) (78) 
i p,{m,}/ 2 ^ , V H 4 V W / / 1 P-{'»/=0} / v ' 

i'G {all legs} v / 

2 ^ W P,{'«/—0}' 1 p,{ra,—0}' ' 

i'G {all legs} 

which holds in the small mass limit up to terms suppressed with the parton masses mf. Of course, 
in the case of massless external lines the respective higher order corrections to the Z-factors in 
Eqs. (1771 ), (|78l ) mostly vanish. Also recall that the amplitude |fW p ) is a vector in color space 
whereas the Z-factors from Eq. © are in this respect simply functions. Non-trivial color depen- 
dence of singularities on the other hand typically arises from soft gluon exchange and therefore 
carries over directly from underlying massless hard scattering amplitude \^ p j mi =o}) ■ Finally, it 
has been emphasized already in the previous discussions, that Eqs. (|76l)-(r78T) require to organize 
the contributions to the massive amplitude \Mpj m .\) in terms of flavor classes, i.e. whether or 

not the heavy parton lines are external. An analogous distinction holds for the gluon factor Z^"' " 1 



when heavy quarks are included for instance as self-energy corrections to the external gluons, see 

,(2) 



Eqs. d89]), (El below. The explicit results for zfV and zf) in Eqs. (g7), (gg) hold for the cases //, 



hi. 

In the light of Eqs. (TT6T) — (TT8T) let us briefly come back to the relation between the factor z( m l°) 
and the perturbative fragmentation functions [19]. Although this connection may come at first as a 
surprise, the two functions are actually intimately related in the context of QCD amplitudes. First 
of all, both are process-independent. Secondly, one may compare both approaches in a computa- 
tion of a one-particle inclusive cross-section of a massive parton based on an amplitude such as 
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Eq. ©. The result takes the form of a convolution of massless cross-section times the perturbative 
fragmentation function. (We refer the reader to the discussion in Refs. [57,68] for complete details 
on this point). Alternatively, we can use Eq. CO) to relate the massive amplitude to the massive 
one. As is clear e.g. from the perturbative expansion in Eqs. (|76|)-(|78T) the proportionality factor 
between \M^s m x) and |^Wp 5 { ra ,=o}) 1S independent of the kinematics and is also not affected by the 
subsequent phase-space integration. Furthermore, this holds separately for virtual and the corre- 
sponding real radiation contributions. Thus, our simple direct relation between massive amplitudes 
and their massless counterpart in Eq. CO) represents the appropriate generalization of the formalism 
of Mele and Nason [19] at the amplitude level. 

In an equivalent formulation, we can also consider the perturbative expansion of Eq. @). To 
that end, we repeat the decomposition of the amplitude |fWpj OT .i) from Eqs. (|76|)-(|78T) up to two 
loops in terms of products of the functions J7 ' p ' , S ' p ' and ' p ' . 



I p,{m,}/ 



\< ] ) 



(i) 

p,{ m i} 

(2) 

P,{ m i} 



> 



) 



I n [ V P (0) >+^ J K j )+K J > 



(0)v 



i'£ {all legs} 



i£ {all legs} 
1 



^-J(W /] ) 2 +^W ]N )Kp (0) > 



(79) 
(80) 

(81) 



■(2)v 



2 . 



€ {all legs} 

where the perturbative expansions of 5^ and |# p ) are defined analogous to Eq. (|75T) . Of course, 
the same qualifications from Section [3] about the distinct flavor classes contributing to the massive 
form factor J also apply here. Now, in the factorization ansatz of Eq. © the function \rt v ) is 
a vector and is a matrix in color space. Thus, their products in Eqs. (l80l) . (I8TI) are in the 
sense of matrix multiplication and all dependence on singular color correlations rests entirely in 
the function S ^ . 

As we remarked above (and as is well known in the literature [11,24,46]) the matrix 5^ is 
subject of a renormalization group equation which allows for an all-order exponentiation of the 
soft contributions. The solution for 5 ^ results in a path-ordered exponential due to mixing of the 
color structures under soft gluon exchange, 



■[p] 



Q 2 

{kij^^si/J 1 ) 



Pexp 



1? 



r^(a(k 2 ,£)) 



(82) 



where P denotes the path ordering. Here, is the so-called soft anomalous dimension, which is 
a matrix in the space of color tensors (see Eqs. ©, ©). Of course, the running coupling a s in the 
argument of is to be taken in d dimensions. For 2 — ► n hard scattering processes with massless 
partons is currently known up to two loops [24,66,69,70] and to one loop results for reactions 
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with massive partons [71,72]. In the latter case, one can show in particular, that the soft anomalous 
dimension r' p ' has a smooth limit for m — ► 0. 

To summarize, we have given in Eqs. (I76l)-(l78l) and (T79l)— <[8TT) two equivalent formulations. 
Both allow to obtain all large logarithms of Sudakov type together with the dimensionally regulated 
soft poles in e and any given QCD amplitude for 2 — > n scattering with parton masses {m,} can 
be constructed by either method. In particular Eq. (l82l) can be used to derive explicit expressions 
for and m Eqs. (|79l)-(l8TI) in terms of the perturbative expansion for the soft anomalous 
dimension r' p l . Most of the other ingredients are explicitly presented in this paper. 

Next, let us discuss the consistency of Eq. © with the results of Ref. [2]. In that reference 
the structure of both soft and collinear singularities for any one-loop amplitude was presented for 
arbitrary values of parton masses. In the approach of Ref. [2] any one-loop rc-parton amplitude can 
be written as: 

|^ p (1) ) = lt\e^ 2 dmj})\^ 0) ) + \^ l)M ), (83) 

where l^p ^) is the Born amplitude for the process under consideration. The amplitude |fM"p^' fin ) 
contains only one-loop corrections which are finite in the limits m ; - — > and £ — > 0. In the following 
we will adapt the results of Ref. [2] to the MS coupling evaluated at a renormalization scale /j. We 
will also assume conventional dimensional regularization for simplicity. In the small mass limit 
the operator In then takes the form (recall that all non-vanishing masses mj are assumed to have 
a common value m): 



J?"W.K}) = t T r T k V Jk (sj k im J ,m*)-£r J ( M ,mfiE) + ...) 



(84) 



where T k are the generators of the gauge group and the kinematical invariants. The dots denote 
mass -independent terms and the functions V ^ are associated to pairs of external partons. One has 
three possible combinations in each pair of partons with two, one or none of them being massive. 
Thus, three separate functions Vji are needed for these three cases. Similarly, the functions Tj are 
different, depending on whether the parton j is massive or massless, i.e. quark, gluon, gluino and 
so on. 

For the sake of comparison with the Z-factor in Eq. ([8]) we write 1/jk and Tj as 

v {2 massive partons) = 2 Ay + (85) 
^(l massive parton) = Ay + ^) (86) 



.(0) 



Vf> = A q + F^ (87) 
rf = n h A g +rf\ (88) 

where the functions Ay, A q and A g are independent of the invariants sjk, i.e. they are the same for 
each external parton (or pair of external partons Therefore one can apply the color algebra to 



3 The superscripts (m) and (0) refer to quantities evaluated in the massive, respectively massless case. 
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express the sum over the products of color generators multiplying these functions directly in terms 
of the corresponding Casimir operators (see Ref. [2]). In this way, all process dependent factors 
are separated into functions that are independent of the mass. All mass dependence on the other 
hand enters only in a process-independent way. Combining the above results one gets 

/i m) (£,A/ 2 ,m 2 ) =/, ( ; 0) ( e y) + £ f q (E,f,m 2 ) + Y,n h f g (E,iiW), (89) 

7=1 7=1 

where /„ is the appropriate operator for purely massless amplitudes [10] evaluated for njr = n { + n h 
light flavors. The function f q is given by one half of the function zj^ (and of course restricted to 
constant terms at order e°) presented in Eq. (1371) . For the function f g we find 

/s (e,^) = -i(i + ln(V)), <*» 

which, when restricted to constant terms at order £°, is related to the function zPj in Eq. (|4Q|) 

through n h f g (£,n 2 ,m 2 ) = zf^/2. 

Finally, we briefly comment on Abelian gauge theories with fermion masses like Quantum 
Electrodynamics (QED). These provide other prominent examples for the considerations of the 
present article. For instance one arrives at QED (with massive electrons) after the usual identifica- 
tion of the color factors, Cp = 1, Ca = and Tf = 1 instead of our QCD convention 7/«y = 
There, the complete calculation of the two-loop radiative photonic corrections in QED to Bhabha 
scattering in the small mass limit have already been performed [20-22]. This included also a 
complete matching at two loops, i.e. the computation of the constant terms which are not logarith- 
mically enhanced. The latter also required the constant terms from the massless, dimensionally 
regularized amplitudes of Ref. [73]. An extension of the results of the present article (and the ex- 
ponentiation in particular) in this direction is a possibility which we leave for a future publication. 



6 Summary 

In this article we have presented a first discussion of the singular behavior of on-shell QCD am- 
plitudes with massive particles beyond one loop. We have performed a systematic study of both, 
the soft singularities typically showing up as poles in £ in dimensional regularization and the struc- 
ture of the large Sudakov (or quasi-collinear) type logarithms of the parton masses, which become 
dominant in the high energy limit. Working in the small mass limit, we have consistently omitted 
power corrections in the parton masses. 

We have presented in Eqs. © and © a general framework for the factorization of n -parton 
amplitudes in QCD which incorporates massive partons. The factorization formula, which we 
have organized in terms of flavor classes, is universal and is valid for any amplitude. We have 
emphasized the strong similarities between scattering amplitudes with massless and massive par- 
tons in the limit where all parton masses are much smaller than the relevant kinematic invariants 
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of the scattering process. In this regime, the factorization formula can be used to directly ob- 
tain (apart from vanishing corrections when the masses tend to zero) the massive amplitude from 
the corresponding massless amplitude, without explicitly computing the former. To that end we 
have introduced the factor z( m '°) as the building block of the proportionality. In the case of heavy 

(m\0) 

quarks we have linked Z}-, to the virtual corrections in the formalism of perturbative fragmen- 
tation function thus generalizing the approach of Ref. [19] to the level of amplitudes. Finally, we 
have explicitly illustrated the predictive power of the factorization ansatz for examples from 2 — > n 
scattering processes in QCD. 

Improved insight into the structure underlying the factorization of amplitudes in the soft and 
(quasi)-collinear momentum regions have enabled us to derive an exponential (1221) for the form 
factor of heavy quarks. We have used this new result to predict the fixed-order expansion of 
the massive form factor to up three loops and, in comparing massless and massive amplitudes, we 
have observed an apparent universality of the respective resummation coefficients G which we find 
worth mentioning. Furthermore, on the basis of Eq. © and the exponentiations for the functions 
3 M and S ' p ' we have shown how to extend our predictions to the perturbative expansion of general 
n-parton amplitudes in QCD with massive partons. 

Thus, the results of the present paper such as Eq. © can be useful to either check explicit eval- 
uations of amplitudes at higher loops or make predictions to higher orders in perturbation theory. 
The material presented can also help to organize calculations, say at NNLO, in terms of divergent, 
but analytically computable, parts and finite remainders that can be integrated numerically. In the 
context of general calculations for differential observables with massive partons at NNLO our fac- 
torization formula may also facilitate the combination of the respective tree-level and one-loop real 
emission amplitudes with the virtual contributions in a process independent manner. 

We will return to these issues as well as potential connections to threshold resummation for 
processes with massive partons in future work. 
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A Appendix 

Here, we give the complete result for one-loop QCD corrections to order £ 2 to the form factor 
of a heavy quark at the scale /u 2 = m 2 in terms of harmonic polylogarithms H mu „^ mw {x), see also 
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Section [3] and Eq. ([431 for definitions. The variable x with < x < 1 for space-like q 2 = —Q 2 < 
is given by 

y/Q2 + 4m 2 -y/& 



a/2 2 + 4m 2 + a/Q 2 " 

For J i we find^ 



(EH) 



5, = ic4-2 + 2(l-^- T l-)H } + C F {-4 + (3- T ^- T ^)Ho 612) 

+ 2(l-^-^)(Ho ) o-2H_ 1) o-C2)}+Bc4-8 + (l- T ^- T ^)(8Ho 
+ 8H i _i 5 o — 4^3 + 4H_i^2 — 4H_i ) o,o — Ho^2 — 4Ho ,-i,o + 2Ho,o,o) 

+ ( 3 " T+i " ill) < H °-° " 2H -''°> " 2 ( 2 " lh ~ TTlM + e2 M- 16 

" 5 ( 4 " T+* " T^)& + ( 3 " T^I " ill) (2H-.C + 4H-,,-,, - 2H_,, ,„ 

— 2^0^,2 — 2Ho ,-i,o+Ho,o,oJ + (l — j — - — ( — ~^2~ + 8H_i^3 — 8H 1 -1^2 

— ^H-^-^-^o + SH-^-^o.o — 16H_i ! o + 2H_i ! o^2 + 8H_i ! o -1,0 — 4H_ 1,0,0,0+ I6H0 
14 

— 3"^ ^ 3 + ^Ho- 1 ^2 + 8 Ho _ 1 _ 1 ,0 — 4Ho _ 1 ,o,o + 8Ho,o — Ho,oC2 — 4Ho,o - 1 ,0 

4 4 



+ 2H„,„,„,„)-2(5- — -— ) fc } 



In Eq. (1AI2I) all harmonic polylogarithms H mu ..^ mw (x), mj = 0, ±1 are understood to be of argu- 
ment x. For the rest, our notation follows Ref. [76] to which the reader is referred for a detailed 
discussion. 

Next we present the one-loop result for the virtual contribution to the perturbative fragmenta- 
tion function to all orders in £ [57]: 

Dr(z) = a s C F ^|±^exp( eTE )r(£) (^V^)- (03) 

As one can easily verify, the expansion of the coefficient of the delta-function in £ coincides 
to all known powers with the factor Zy? in Eq. (|37|) . which suggests that Eq. (1AI3I) is indeed the 
proper generalization to all orders in £. 

Finally, we discuss the derivation of the one-loop heavy quark insertion in the tree-level gluon 
form factor. It is clear that this diagram coincides with the one loop contribution to . We find 
it particularly convenient to evaluate this diagram in a physical light-cone gauge. Following the 
procedure outlined in Refs. [57, 58] we obtain to all orders in £: 

Zg } = a s n h (-^) exp(£YE)r(£) (^) . ©4) 



Upon expansion in £ Eq. (1401) is derived. 

4 We thank J. Gluza for providing us with the integral SE212m of Refs. [12, 13] to order e 3 , see also 

|http : / /www-zeuthen .desy . de/theory/ research/bhabha/ | 
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